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Abstract 

We  consider  the  satisfiability  problem  for  Boolean  combinations  of  generalized  2SAT  constraints,  which  are 
linear  constraints  with  at  most  two,  possibly  unbounded,  integer  variables  having  coefficients  in  {—1, 1}. 
We  prove  that  if  a  satisfying  solution  exists,  then  there  is  a  solution  with  each  variable  taking  values  in 
[— n  •  (6max  +1)5  n  ■  (6max  +  1)]>  where  n  is  the  number  of  variables,  and  6max  is  the  maximum  over  the 
absolute  values  of  constants  appearing  in  the  constraints.  This  solution  bound  improves  over  previously 
obtained  bounds  by  an  exponential  factor.  Our  result  enables  a  new  enumerative  approach  to  satisfiability 
checking.  An  experimental  evaluation  demonstrates  the  efficiency  of  this  approach  over  previous  techniques. 
As  a  corollary  of  our  main  result,  we  obtain  a  polynomial-time  algorithm  for  approximating  optima  of 
generalized  2SAT  integer  programs  to  within  an  additive  factor. 
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1  Introduction 


Generalized  2SAT  constraints  are  a  special  class  of  linear  constraints  over  integer  variables.  A  generalized 
2SAT  (G2SAT)  constraint  (also  called  a  unit  two  variable  per  inequality  or  UTVPI  constraint)  has  at  most 
two  variables,  and  variable  coefficients  are  in  {—1, 1}.  The  variables  are  not  required  to  have  finite  upper 
or  lower  bounds.  Useful  optimization  problems,  such  as  the  minimum  vertex  cover  and  the  maximum 
independent  set  problems,  can  be  modeled  using  generalized  2SAT  constraints,  and  several  applications 
of  constraint  logic  programming  and  automated  theorem  proving  also  generate  G2SAT  constraints  (e.g., 
see  [14,  1]). 

A  G2SAT  formula  is  a  Boolean  combination  of  G2SAT  constraints.  In  this  paper,  we  consider  the 
problem  of  checking  the  satisfiability  of  G2SAT  formulas.  It  is  easily  seen  that  this  problem  is  NP-complete. 
However,  the  special  case  of  checking  satisfiability  of  a  conjunction  of  G2SAT  constraints  (i.e.,  finding  a 
feasible  integer  poinf  in  a  G2SAT  polyhedron)  can  be  solved  in  polynomial  lime;  for  example,  a  modified 
version  of  Fourier-Mofzkin  eliminalion  [8,  22]  (reviewed  in  Seclion  2)  runs  in  O(n^)  lime. 

Currenl  approaches  (e.g.,  [1])  lo  checking  Ihe  salisfiabilily  of  a  G2SAT  formula  employ  a  combination 
of  Boolean  salisfiabilily  solving  and  linear  conslrainl  solving.  Trulh  values  are  assigned  fo  linear  conslrainls 
so  lhal  Ihe  G2SAT  formula  is  satisfied.  Each  such  Irulh  assignmenl  corresponds  lo  a  G2SAT  polyhedron.  If 
Ihis  polyhedron  has  a  feasible  integer  poinf,  lhal  poinf  satisfies  Ihe  original  G2SAT  formula  as  well.  If  nol, 
anolher  Irulh  assignmenl  musl  be  found.  Given  a  G2SAT  formula  f  wilh  m  conslrainls  and  n  variables,  and 
assuming  lhal  integer  feasibility  is  checked  using  Ihe  afore-menlioned  modified  Fourier-Mofzkin  elimination 
algorilhm,  Ihe  currenl  techniques  have  a  worsl-case  running  time  of  0(2"*  •  n^).^ 

In  Ihis  paper,  we  prove  lhal  a  satisfying  solution  exisls  for  a  G2SAT  formula  f  if  and  only  if  Ihere  is  a 
solution  lo  (j)  wilh  each  variable  faking  values  in  Ihe  finite  range  [— n  •  (6max  +  I),”  •  (&max  +  1)]>  where 
n  is  Ihe  number  of  variables  in  f,  and  bmax  is  the  maximum  over  the  absolute  values  of  constant  terms  in 
the  constraints.  That  such  a  bounded  solution  exists  is  not  surprising,  since  satisfiability  solving  of  G2SAT 
formulas  is  in  NP.  However,  the  previously  best  known  solution  bounds  [4,  23,  15,  18]  are  H(n^  •  (6max  + 
1)  •  2").  In  particular,  our  result  eliminates  the  2"  term,  thereby  exponentially  reducing  the  solution  bound. 

Our  result  can  be  used  to  check  satisfiability  of  G2SAT  formulas  in  worst-case  time  0(2"'°s<t)  .^^here 
d  =  2  ■  n  ■  (bmax  +  1)>  by  encoding  each  integer  variable  with  logd  Boolean  variables.  This  yields  a  more 
efficient  satisfiability  checker  for  highly  over-constrained  formulas,  where  m  =  H(n  •  logd).^  The  latter  is 
often  the  case  for  theorem  proving  applications  in  program  analysis  and  hardware  verification. 

A  key  step  in  our  proof  is  to  show  that  for  a  G2S  AT  polyhedron,  if  a  feasible  integer  point  exists,  then  one 
exists  within  a  unit  hypercube  centered  at  any  minimal  face  solution  (extreme  point).  As  a  corollary  of  this 
result,  we  obtain  a  polynomial-time  algorithm  for  approximating  optima  to  an  additive  factor  in  generalized 
2SAT  integer  programs. 

Our  theoretical  results  are  validated  by  an  experimental  evaluation  on  randomly  generated  G2SAT  for¬ 
mulas,  which  shows  that  a  decision  procedure  based  on  our  approach  can  greatly  outperform  other  proce¬ 
dures. 

1.1  Related  Work 

There  has  been  much  previous  work  on  integer  programming  with  two  variables  per  inequality  (see,  e.g., 
the  work  by  Hochbaum  et  al.  [13,  12,  11]).  The  main  differences  between  this  work  (applied  to  G2SAT 
constraints)  and  ours  are  threefold.  First,  our  focus  is  on  satisfiability  solving  of  arbitrary  G2SAT  formulas 

'Assuming  the  trivial  worst-case  bound  of  0{2^)  for  checking  satisfiability  of  a  Boolean  formula  in  N  variables. 

^For  a  conjunction  of  G2SAT  constraints,  m  is  O(n^),  since  one  can  eliminate  redundant  constraints.  However,  for  an  arbitrary 
Boolean  combination,  this  is  not  the  case. 
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and  not  linear  optimization  over  G2SAT  polyhedra.  Second,  we  do  not  require  variables  to  be  bounded. 
Finally,  for  our  approximation  result,  the  objective  function  can  be  an  arbitrary  linear  function,  without  any 
restriction  on  the  sign  of  cost  coefficients. 

Previous  results  on  bounding  solutions  have  been  derived  in  the  context  of  showing  that  integer  linear 
programming  is  in  NP  [4,  23,  15,  18].  Even  when  specialized  for  G2SAT  integer  programs,  these  bounds 
are  •  (6max  +  1)  •  2”) .  Our  result  is  therefore  an  exponential  reduction  in  the  solution  bound  for  G2S AT 
integer  programs,  and,  to  the  best  of  our  knowledge,  has  not  been  obtained  before. 

Our  results  rely  on  the  modified  version  of  Fourier-Mofzkin  elimination  for  checking  infeger  feasibilily 
of  a  G2SAT  polyhedron;  fhis  algorifhm  is  described  by  Subramani  [22],  and  an  incremenfal  version  has  been 
given  by  Harvey  and  Sfuckey  [10]. 

Theorem  provers  fhaf  can  check  G2SAT  formulas,  such  as  CVC-Life  [7],  are  essentially  a  combination 
of  a  SAT  solver  and  a  solver  for  a  system  of  linear  consfrainfs.  In  fhe  case  of  CVC-Life,  fhis  solver  is 
fhe  Omega  fesf  [19],  which  for  G2SAT  consfrainfs  is  idenfical  fo  fhe  modified  Fourier-Mofzkin  elimination 
algorithm  referenced  above. 

1.2  Outline 

The  rest  of  the  paper  is  organized  as  follows.  We  begin  in  Section  2  with  useful  background  definitions  and 
results.  Section  3  contains  the  main  theoretical  contributions  of  this  paper.  We  present  experimental  results 
in  Section  4  and  conclude  in  Section  5. 

2  Background 

We  state  here,  in  brief,  some  definitions  and  theorems  used  in  the  remainder  of  the  paper.  Further  details 
can  be  found  in  standard  textbooks  on  polyhedral  theory  and  integer  linear  programming  (e.g.,  [17,  20]). 

Following  standard  linear  programming  notation,  we  denote  the  number  of  variables  by  n  and  number 
of  constraints  by  m.  We  assume  that  a  linear  constraint  is  specified  in  fhe  form  a  •  x  >  b,  where  a  is 
a  n-dimensional  infeger  vector  [oi,  02,  •  •  •  ,  On],  x  is  a  n-dimensional  vecfor  of  infeger-valued  variables 
[xi,X2,  -  ■  ■  ,  Xn],  and  b  is  an  infeger.  A  sysfem  of  consfrainfs  is  specified  as  A  •  x  >  b,  where  A  is  a  m  x  n 
mafrix  wifh  integral  enfries,  b  is  a  m  x  1  integer  vector  [bi,b2,  ■  ■  ■  ,  and  x  is  a  n  x  1  vecfor  of  infeger- 
valued  variables.  We  use  bmax  to  denote  fhe  Lao  norm  of  b;  i.e.,  bmax  =  max,  |6j|. 

The  terms  feasible  and  satisfiable  are  used  interchangeably,  as  also  are  lattice  point  and  integer  point. 

2.1  G2SAT  Formulas 

Definition  2.1  A  constraint  a  -  >  b  is  said  to  be  an  absolute  constraint  if  exactly 

zero,  a  difference  constraint  if  exactly  two  of  the  OiS  are  non-zero  with  one  being  -|-1 
a  sum  constraint  if  exactly  two  of  the  aiS  are  non-zero  with  both  -|-1  or  both  —1. 
a  -  5t>b  is  said  to  be  a  G2SAT  consfrainf  if  it  is  either  an  absolute,  a  difference, 

A  G2SAT  formula  is  generafed  by  fhe  following  grammar: 

(j)  ::=  true  |  false  \xi X2>b\xi  —  X2>b\x  >b  \  |  A  ^2  |  V  ^2 

Notice  fhaf  a  negation  on  a  G2SAT  consfrainf  can  be  eliminafed  by  rewriting  fhe  consfrainf.  A  G2SAT 
consfrainf  remains  G2SAT  under  such  rewrifing.  The  only  change  is  to  fhe  sign  of  variable  coefficienls,  and 
to  fhe  consfanf  ferm,  which  can  increase  in  absolute  value  by  af  mosf  1. 


one  of  the  aiS  is  non- 
and  the  other  —  1,  and 

or  a  sum  constraint. 
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Example  2.1  Consider  the  following  G2SAT  formula 

(-1  xi-\-  X2>  —1)  A  {x2  —  >  0  V  Xi>\^ 

The  constraint  xi  +  X2  >  is  a  sum  constraint,  X2  —  X3  >  0  is  a  difference  constraint,  and  3^4  >  1  is  an 
absolute  constraint.  The  negation  can  be  eliminated  to  obtain  an  equivalent  G2  SAT  formula 

—xi  —  X2  >2  A  [x2  —  X3  >  0  V  3:4  >  1) 

Note  that  the  value  o/6max  htis  increased  from  1  to  2  after  eliminating  the  negation. 

Not  all  families  of  linear  constraints  are  closed  under  eliminating  negations.  For  example,  the  class  of 
Horn-SAT  constraints,  which  comprises  all  constraints  with  at  most  one  variable  with  a  positive  coefficient, 
are  not  closed  under  eliminating  negations. 

Definition  2.2  Given  a  G2  SAT  formula  f,  an  enumeration  bound  is  an  integer  d  such  that  f  is  lattice  point 
feasible  if  and  only  if  it  contains  a  lattice  point  in  the  n-dimensional  hypercube  11^=1  interval 

[— d,  d\  is  termed  as  an  enumeration  domain. 

2.2  Polyhedral  Theory 

Definition  2.3  A  minimal  face  of  a  polyhedron  is  a  face  that  does  not  contain  any  other  face  of  the  polyhe¬ 
dron.  A  point  lying  on  a  minimal  face  is  called  a  minimal  face  solution  (MFS). 

When  the  minimal  face  is  an  extreme  point,  a  MFS  is  a  basic  feasible  solution. 

We  write  (A',  b')  C  (A,  b)  to  indicate  that  the  polyhedral  system  A'  •  x  >  b'  is  a  subsystem  of  the 
polyhedral  system  A  •  x  >  b.  Also,  for  a  matrix  A,  let  r(A)  denote  the  rank  of  A.  We  have  the  following 
characterization  of  a  minimal  face. 

Theorem  2.1  ([20])  Let  P  =  {x  :  A  •  x  >  b}  denote  a  polyhedron.  A  non-empty  subset  E  CD  is  a  mini¬ 
mal  face  ofP,  if  and  only  ifF  =  {x  :  A'  •  x  =  h'},  for  some  system  A'  •  x  >  b',  where  (A',  b')  C  (A,  b), 
andr{A'^  b')  =  r(A,  b). 

Fourier-Motzkin  (FM)  elimination  [8]  is  a  well-known  projection  technique  for  polyhedra.  Starting  with 
a  polyhedron  P  :  A  •  x  >  b,  a  variable  Xj  is  projected  using  FM  elimination  in  the  following  steps: 

1.  Partition  the  system  of  constraints  into  three  sets  Pj,  Nj,  Zj  as  follows.  For  each  constraint  i,  1  < 
i  <  m,'we  add  it  to: 

Pj,  if  Q.i,j  ^  O9 
Nj,  if  Oij  <  0; 

Zj,  otherwise. 

2.  Initialize  the  set  of  new  constraints,  $,  to  Zj. 

3.  For  every  pair  of  constraints  (ip,  ip/),  where  ip  E  Pj  and  ip/  E  Nj,  add  the  following  constraint  to  $: 

n 

{fiip,3  ■  ^iN,k  ~  Clip,k  •  CliN,j)  '  ^  h 

k=l 
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Clearly,  the  eoeffieient  of  Xj  in  every  eonstraint  in  $  is  0.  Denote  the  polyhedron  defined  by  the  set  $  as 

P  :  A  •  X  >  b.  It  is  easy  to  see  that  P  has  a  solution  in  M”  if  and  only  if  P  has  a  solution  in  . 

Note  that  if  P  is  G2S  AT,  P  need  not  be  G2SAT.  However,  it  is  possible  to  modify  the  basie  FM  elimina¬ 
tion  proeedure  by  adding  a  coefficient  normalization  step,  so  that  the  resulting  polyhedron  remains  G2SAT, 
and  moreover,  is  lattiee  point  feasible  iff  P  is.  Notiee  that  the  only  non-G2SAT  eonstraints  in  P  are  of  the 
form  2xi  >  6  or  —2xi  >  b,  obtained  by  adding  a  sum  eonstraint  involving  Xi  and  Xj  with  a  differenee 
eonstraint  involving  those  variables.  By  dividing  both  sides  of  a  newly  ereated  non-G2SAT  eonstraint  by  2, 
and  rounding  up  the  RHS  if  it  is  an  odd  multiple  of  we  obtain  a  G2SAT  eonstraint  with  the  same  integral 
solutions  as  the  original.  In  this  way,  we  replaee  eaeh  non-G2SAT  eonstraint  in  P  with  a  eorresponding 
G2SAT  eonstraint  to  obtain  a  G2SAT  polyhedron  P'  :  A'  •  x'  >  b'. 

We  will  refer  to  the  modified  FM  elimination  proeedure  as  Fourier-Motzkin  elimination  with  coefficient 
normalization  (FM-CN).  If  is  easy  fo  see  fhaf  FM-CN  preserves  infegral  solutions,  i.e.,  P  is  laffiee  poinf 
feasible  iff  P'  is.  One  ean  use  FM-CN  fo  eheek  fhe  feasibilify  of  G2SAT  polyhedra  in  time  O(n^),  by 
sueeessively  eliminafing  variables,  eheeking  af  eaeh  sfep  fhaf  we  do  nof  generafe  a  frivially  false  eonsfrainf. 
Af  any  sfep,  we  are  guaranfeed  fo  have  a  sysfem  of  no  more  fhan  0{v?)  eonsfrainfs,  sinee  fhere  are  only 
4  •  (2)  possible  non-redundanf  G2SAT  eonsfrainfs  on  n  variables. 

3  Theoretical  Results 

Our  fheorefieal  resulfs  are  organized  as  follows.  We  begin,  in  Seefion  3.1,  by  showing  fhaf  if  a  G2SAT 
polyhedron  has  a  minimal  faee  solution  (MFS),  fhen  fhere  exisfs  a  MFS  wifh  eaeh  eomponenf  half-infegral 
and  in  [— n  •  femax,  n  ■  &max]-  The  main  fheorem,  presenfed  in  Seefion  3.3,  enables  us  fo  go  from  bounding  a 
MFS  fo  bounding  infeger  solufions.  This  fheorem  sfafes  fhaf  if  a  G2SAT  polyhedron  is  infeger  feasible,  fhen 
if  is  possible  fo  find  a  infegral  solution  wifhin  a  unif  box  eenfered  af  any  MFS;  i.e.,  by  “rounding”  a  MFS. 
In  fhis  seefion,  we  also  deseribe  how  fo  exfend  resulfs  for  G2SAT  polyhedra  fo  arbifrary  G2SAT  formulas. 
Seefion  3.2  presenfs  auxiliary  resulfs  on  rounding  fhaf  are  used  fo  prove  fhe  main  fheorem.  Finally,  in 
Seefion  3.4,  we  show  fhaf  fhe  main  fheorem  ean  be  used  fo  obfain  an  additive  approximafion  resulf  for 
opfimizing  an  arbifrary  linear  eonsfrainf  over  a  G2SAT  polyhedron. 

3.1  Minimal  Face  Solutions  of  G2SAT  Polyhedra 

We  begin  wifh  a  useful  lemma. 

Lemma  3.1  Let  P  :  A  •  x  >  b  represent  a  system  ofm  difference  constraints  on  n  variables.  Then,  P  has 
a  feasible  integer  solution  if  and  only  if  it  has  an  integer  solution  in  the  hypercube  nr=i[o,  (n-i)-  ^max]- 

Proof:  Consider  fhe  dual  eonsfrainf  graph  as  ouflined  by  Cormen  ef  al.  [5].  A  solution  fo  fhe  sysfem  is 
obfained  by  assigning  fo  eaeh  variable,  fhe  shorfesf  pafh  from  fhe  souree.  The  lengfh  of  any  shorfesf  pafh  is 
bounded  by  (n  —  1)  •  6max-  The  resulf  follows.  □ 

The  following  lemma  eonsiders  bounding  a  MFS  of  a  G2SAT  polyhedron  in  fhe  non-negafive  orfhanf. 

Lemma  3.2  Let  P:A-x>b,  x>0  denote  an  arbitrary  G2SAT  polyhedron  in  the  non-negative  orthant 
with  m  constraints  and  n  variables.  Then,  if  a  MFS  exists,  then  there  is  a  MFS  with  each  component 
half-integral  and  at  most  n  ■  6max- 

Proof:  Suppose  polyhedron  P  has  a  minimal  faee  solution.  Hoehbaum  ef  al.  [13]  have  shown  fhaf  fhis  MFS 
musf  be  half-infegral.  We  foeus  here  on  showing  fhe  n  •  6max  bound. 
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By  definition,  the  minimal  faee  eorresponding  to  this  MFS  satisfies  a  system  A'  •  x  =  b',  where 
(A'  b')  C  (A  b),  and  r{A')  =  r(A)  =  k  for  some  1  <  A:  <  n  (assuming,  w.l.o.g.,  that  m  <  n).  Ae- 
eordingly,  there  are  k  independent  variables  and  n  —  k  dependent  variables  in  the  system;  without  loss  of 
generality,  we  assume  that  the  first  k  variables  are  independent  and  set  the  dependent  variables  to  0.  This 
results  in  a  system  Pi  :  A"  •  x"  =  b",  x"  >  0,  where  the  eomponents  of  b"  are  also  eomponents  of  b,  and 

X"  =  [xi,X2,... 

The  system  Pi  eontains  3  types  of  eonstraints  (equations),  viz.,  absolute,  differenee  and  sum.  We 
eonsider  eaeh  of  these  types  in  turn: 

1.  An  absolute  eonstraint  is  of  the  form  Xi  =  b.  Sinee  x"  >  0,  the  value  of  Xi  must  be  in  [0,  6max]- 

2.  A  sum  eonstraint  ean  be  written  in  the  form  Xi  +  Xj  =  b,  where  6  >  0.  Sinee  x"  >  0,  it  follows  that 

0  ^  5  Xj  ^  b  ^  bjnax  ■ 

3.  From  the  two  eases  above,  we  eonelude  that  the  value  of  any  variable  appearing  in  an  absolute  or  sum 
eonstraint  must  lie  in  [0, 6max]  (and  moreover,  there  exists  sueh  a  half-integral  value). 

W.l.o.g,  let  xi,  X2,  ■  ■  ■  ,  Xi,  I  <  k,be  variables  appearing  in  the  absolute  and  sum  eonstraints,  and 
let  x\,X2t  ■  ■  ,x’l  be  the  eorresponding  half-integral  values  in  [0,  femax]  satisfying  these  eonstraints. 
Substituting  these  values  into  the  differenee  eonstraints  might  ereate  new  absolute  eonstraints,  but  no 
new  differenee  or  sum  eonstraints.  The  eonstant  term  in  new  absolute  eonstraints  generated  thus  is 
half-integral  and  of  absolute  value  at  most  26niax-  The  substitution  proeess  ean  be  iterated  at  most 
A:  —  1  times  leading  to  absolute  eonstraints  with  half-integral  eonstant  terms  at  most  k  ■  6max-  Thus, 
a  variable  appearing  in  any  of  the  absolute  eonstraints  generated  in  this  iterative  proeess  takes  half¬ 
integral  values  in  [0,  A:  •  6max]- 

When  the  above  iterative  substitution  proeess  terminates,  the  only  eonstraints  possibly  left  are  some 
of  the  original  differenee  eonstraints,  eaeh  with  an  integral  eonstant  term  of  absolute  value  at  most 
bjnax-  Sinee  these  eonstraints  are  satisfiable,  we  ean  apply  Lemma  3.1  to  eonelude  that  there  exists  a 
solution  to  these  eonstraints  with  eaeh  variable  taking  integral  values  in  [0,  (A:  —  1)  •  6max]  (sinee  at 
most  k  variables  appear  in  these  eonstraints). 

Sinee  k  <  n,  we  eonelude  that  there  exists  a  solution  to  Pi  with  eaeh  eomponent  at  most  n  •  6max-  LI 

We  now  generalize  the  result  to  an  arbitrary  G2SAT  polyhedron. 

Theorem  3.1  LetF  :  A  ■  X  >  b  denote  an  arbitrary  G2SAT polyhedron  with  m  constraints  and  n  variables. 
If  a  MFS  exists,  there  exists  a  MFS  with  each  component  half-integral  and  in  the  interval  [-n-bmax,  n'bmax]- 

Proof:  Suppose  X*  is  a  MFS  of  P.  Letji,j2,...  ,  be  the  set  of  all  eolumn  indiees,  1  <  ^jk  ^  ft, 

sueh  that  x*j^  <  0  for  all/,  1  <  t  <  k.  Construet  a  matrix  A'  by  multiplying  the  j)th  eolumn  of  A  by  —1  for 
all  I,  leaving  other  eolumns  unehanged.  We  observe  that: 

1.  The  polyhedron  P'  :  A'  •  x  >  b,  x  >  0  is  also  G2SAT. 

2.  If  we  eonstruet  x'*  from  x*  by  negating  x*^  for  all  I,  1  <  I  <  k,  x.'*  satisfies  P'.  Moreover,  we  argue  that 
it  is  a  MFS  of  P'  as  follows: 

Let  (A,b)  C  (A,b)  be  the  eonstraints  satisfied  with  equality  at  x*,  and  (A',b')  C  (A',b')  be  the 
eonstraints  satisfied  with  equality  at  x'*.  Then,  r(A)  =  r(A'),  sinee  A  and  A'  eorrespond  to  the  same 

rows  (of  A  and  A'  respeetively).  Also,  note  that  r(A)  =  r(A').  Finally,  sinee  (A,b)  define  a  minimal 
faee  of  P,  r(A)  =  r(A)  [20]. 

Thus,  r(A')  =  r(A'),  and  so  x'*  is  a  MFS  of  P'. 
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Using  an  identical  argument,  we  conclude  that,  from  a  MFS  of  P',  we  can  construct  a  MFS  of  P  by  negating 
values  to  Xj^ ,  Xj^ ,  ■  •  •  ,  Xj^, . 

Since  P'  has  a  MFS,  by  Lemma  (3.2)  it  must  have  a  MFS  with  each  component  half-integral  and  in 
[0,n-6max]-  It  follows  that  P  has  a  MFS  with  each  component  half-integral  and  in  [— n  •  femax,  n  ■  6max]- 


3.2  Rounding  and  Semi-Rounding 

Definition  3.1  A  rational  number  x  is  said  to  be  odd  half-integral  if  it  is  an  odd  multiple  of 

Definition  3.2  A  vector  z  is  said  to  be  a  rounding  of  a  vector  Stifz  is  integral  and  |  |z  —  x|  |oo  <  j- 

Definition  3.3  A  vector  z  is  said  to  be  a  semi -rounding  of  a  vector  x  if  all  of  the  following  conditions  hold: 
(1)  ||z  —  x||oo  <  5/  (2)  all  components  ofz  are  half-integral;  and  (3)  if  a  component  of  5l  is  integral,  so  is 
the  corresponding  component  of  z. 

Lemma  3.3  Let  a  -  x  >b  be  a  G2SAT  constraint.  Let  x*  be  a  half-integral  vector  such  that  a  •  x*  >  b,  and 
let  w*  be  an  arbitrary  semi-rounding  ofx*.  Then,  a  •  w*  >  b. 

Proof:  The  proof  proceeds  by  case  splitting  on  the  number  of  variables  in  the  constraint. 

1.  Suppose  the  constraint  involves  only  one  variable.  Then,  it  is  either  of  the  form  Xi  >  b  or  —Xi  >  b. 

Correspondingly,  we  either  have  x^  >  b  or  — >  b.  Since  x^  is  half-integral,  in  both  cases  the  LHS 
exceeds  b  by  at  least  Thus,  any  semi -rounding  of  x^  satisfies  the  constraint. 

2.  Suppose  the  constraint  has  two  variables,  Xi  and  Xj.  Then,  since  x^  and  x*  are  both  half-integral,  one 
of  the  following  two  cases  must  hold: 

(a)  The  LHS  is  integral,  and  exceeds  b  by  at  least  1.  But  any  semi-rounding  of  x^  and  x*  can 
decrease  the  LHS  by  at  most  1,  and  hence  satisfies  fhe  consfraint. 

(b)  The  LHS  is  odd  half-integral,  i.e.,  one  of  x^  and  Xj  is  integral  and  the  other  odd  half-integral. 

Thus,  the  LHS  exceeds  6  by  at  least  In  this  case,  any  semi-rounding  of  x^  and  x*  can  decrease 

the  LHS  by  at  most  j,  and  will  satisfy  the  constraint. 

□ 

Since  every  rounding  z  of  x*  is  also  a  semi-rounding  of  x*,  we  obtain  the  following  corollary: 

Corollary  3.1  Let  a  -  5i  >  b  be  a  G2SAT  constraint.  Let  x*  be  a  half-integral  vector  such  that  a  •  x*  >  b, 
and  let  z  be  an  arbitrary  rounding  ofx*.  Then,  a  -  z  >  b. 

We  now  state  a  useful  property  of  Fourier-Motzkin  elimination  with  coefficient  normalization. 

Proposition  3.1  Let  P  :  A  •  x  >  b  denote  a  G2SAT  polyhedron  in  and  x*  = 

denote  a  half-integral  feasible  solution  to  P.  Further,  suppose  that  P  is  lattice  point  feasible. 

Let  P'  :  A'  •  x'  >  b'  be  obtained  from  P  by  projecting  out  variable  Xn+i  using  Fourier-Motzkin  elimi¬ 
nation  with  coefficient  normalization  and  denote  {x\^X2t  ■  ■  ^xf)  byx!* .  Then,  there  exists  a  semi-rounding 
w'*  o/x'*  such  that  w'*  is  a  solution  to  P^ 

Proof:  First,  note  that  since  P  is  lattice  point  feasible,  so  is  P'. 

If  x'*  is  already  a  solution  to  P'  then  the  theorem  holds  trivially. 

So  suppose  that  x'*  does  not  satisfy  P'.  The  only  reason  this  occurs  is  because  x'*  is  cut  off  by  coefficient 
normalization,  i.e.,  due  to  the  presence  of  one  or  both  of  the  following  situations: 
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1.  There  exists  at  least  one  variable  Xi,i  E  I,  such  that  P  has  constraints  of  the  form: 


(1) 

(2) 


Xi  +  Xn+1  >  h\ 


which  result  in  the  following  constraint  in  P': 


Xi  > 


'bi  +  K 
2 


where,  6,  +  6'  is  odd. 

Since  x'*  does  not  satisfy  P',  the  following  equality  also  holds: 

x*  =  h±A 

i  o 


2.  There  exists  at  least  one  variable  xj,  j  E  J,  sueh  that  P  has  eonstraints  of  the  form: 


-Xj  -  Xn+1  >  b'j 


which  result  in  the  following  constraint  in  P': 


Xj  < 


where,  bj  +  b'j  is  odd. 

Since  x'*  does  not  satisfy  P',  the  following  equality  also  holds: 


(3) 


(4) 

(5) 

(6) 

(V) 


(8) 


Note  that  for  some  i  E  I,  and  j  E  J,iii  =  j,  then  we  must  have  But  that  would  mean  that 

P'  is  infeasible,  since  constraints  (3)  and  (7)  would  contradict  each  other.  Hence,  we  can  assume  hereafter 
that  the  two  index  sets  I  and  J  are  disjoint. 

We  now  give  a  rounding  algorithm  that  generates  a  semi-rounding  w'*  of  x'*  that  satisfies  P'.  The 
rounding  algorithm  is  as  follows: 

1.  Initialize  the  set  of  variables  to  be  rounded  up,  14,  to  be  {xi\i  E  I}.  Similarly,  initialize  the  set  of 
variables  to  be  rounded  down,  T>  as  {xj\j  E  J}. 

2.  Uo  ■=  U,  Vo  :=  V,  t  :=  0. 

3.  Compute  Ut+i  and  Vt+i  as  follows.  For  every  Xi  E  lAt  and  Xj  E  Vt, 

(a)  Include  in  Ut+i  any  variable  Xk  such  that  the  following  constraints  in  P',  which  are  valid  for  P, 
hold  with  equality  at  x'* : 


^  bfi^i 

4”  ^k  —  bjk 
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(9) 

(10) 


(b)  Include  in  Vt+i  any  variable  Xk  such  that  the  following  constraints  in  P',  which  are  valid  for  P, 
hold  with  equality  at  x'* : 


-Xk-Xi  >  b'l^i  (11) 

Xj  -  Xk  >  (12) 

4.  If  Ut+i  C  U  and  22t+i  C  V,  stop. 

Otherwise,  perform  the  assignments  U  :=UU  lAt+i,  V  :=  T>U  Vt+i,  f  :=  f  +  1,  and  go  to  step  (3). 

It  is  easy  to  prove  by  induction  on  t,  that  for  any  Xk  E  U,  k  ^  I,  there  either  exists  i  E  I  and  an  integer  bki 
such  that 


4  -  4  =  hi  (13) 

or  a  j  E  J  and  an  integer  bjk  such  that 

x*  +  x*k  =  bjk  (14) 

Similarly,  for  each  Xk  E  V,  k  ^  J,  there  either  exists  i  E  I  and  an  integer  b'l^-  such  that 

-x*k  -  4  =  b'ki  (15) 

or  a  j  E  J  and  an  integer  such  that 

x*-xl  =  b'jk  (16) 

Suppose  the  two  sets  lA  and  V  are  disjoint.  Then,  to  obtain  a  semi-rounding  w'*  of  x'*,  we  round  up 
every  variable  in  U  and  round  down  every  variable  in  T). 

To  complete  the  proof,  the  following  two  sub-goals  remain  to  be  established: 

1.  U^V  =  %. 

2.  w'*  satisfies  P'. 

Assuming  the  first  sub-goal,  consider  the  second  sub-goal  first.  We  observe  that: 

•  By  Lemma  3.3,  any  constraints  in  P'  that  are  not  satisfied  with  equality  at  x'*  will  continue  to  be 
satisfied  by  w'*. 

•  From  Equations  (13)-(16),  we  nofe  that  for  all  Xk  E  U  U  V,  x*f,  is  odd  half-integral,  since  it  is  an 
integral  offset  from  x^  or  Xj  for  some  i  E  I  or  j  E  J. 

Thus,  for  all  Xk  EUUV,  there  cannot  be  any  absolute  constraint  involving  Xk  in  P'  that  holds  with 
equality  at  x'*.  Thus,  by  Lemma  3.3,  the  semi-rounding  produced  by  the  above  algorithm  satisfies 
these  absolute  constraints. 

•  Steps  3(a)  and  3(b)  of  the  rounding  algorithm  ensure  that  all  two-variable  constraints  of  P'  satisfied 
with  equality  at  x'*  continue  to  be  satisfied  by  the  generated  semi-rounding.  For  example,  if  Xk  —  Xi  > 
bki  is  satisfied  with  equality  at  x'*,  and  x^  is  rounded  up,  so  is  x*^,,  so  the  constraint  continues  to  be 
satisfied. 


Thus,  if  the  two  sets  U  and  V  are  disjoint,  we  ean  eonelude  that  w'*  satisfies  P'.  We  will  now  show  that 
the  former  is  indeed  the  ease. 

The  proof  is  by  eontradietion.  Suppose  W  fl  7^  0.  Let  Xk  be  a  variable  present  in  both  sets.  As  we 
noted  before,  for  any  ?  E  /  and  j  E  J,  i  ^  j,  so  we  ean  assume  that  k  is  neither  in  I  nor  in  J.  We  have  the 
following  eases,  eaeh  of  whieh  leads  to  a  eontradietion: 

1.  Equations  (13)  and  (16)  hold.  Then,  for  some  integer  bji,  we  have 

X*  -  X*  =  bji  (17) 

The  above  equation  eorresponds  to  the  following  inequality  derived  by  adding  Inequalities  (9)  and  (12), 
whieh  is  valid  for  both  P  and  P': 


Xi  ^  bji 


(18) 


Further,  from  Equation  (17)  and  Inequalities  (1), 

-bji  =  X*  - 
-bji  =  X*  - 

Also  from  Equations  (4)  and  (8),  we  know  that 


From  (19),  (20),  and  (21)  above,  we  infer  that  bi  +  bj  =  &^  +  b'j  =  —bji. 

Thus,  the  inequalities  in  (19)  and  (20)  hold  with  equality.  Also,  from  Inequalities  (1)  and  (5),  Xi—Xj  > 
bi  +  bj  is  valid  for  P.  Thus,  we  ean  eonelude  that  Inequality  (18)  holds  with  equality  for  P.  This 
further  implies  that  Inequalities  (1),  (2),  (5),  and  (6)  hold  with  equality  for  P. 

Sinee  there  is  a  unique  solution  to  Constraints  (1),  (2),  (5),  (6)  and  (18)  that  satisfies  them  with  equal¬ 
ity,  in  every  feasible  solution  of  P,  Xj  =  x^,  Xj  =  Xj,  and  Xn+i  =  2:^-1- 1-  Sinee  at  least  one  of  x^  and 
X*  is  odd  half-integral,  this  eontradiets  the  premise  that  P  has  a  lattiee  point  solution. 


(2),  (5),  and  (6),  we  ean  eonelude  that 


>  +  bj 

(19) 

IV 

+ 

(20) 

+  6, +  &'  +  (/, 

2.  Equations  (14)  and  (15)  hold.  This  ease  is  identieal  to  Case  (1)  above. 

3.  Equations  (14)  and  (16)  hold.  Then,  we  have 


bjk  +  ^jk 
2 


(22) 


bq  1,  —  bo  —  b^- 

This  implies  that  — . 

Further,  Equation  (22)  eorresponds  to  the  following  valid  eut  for  P'  (i.e.,  it  preserves  lattiee  point 
solutions),  obtained  by  adding  (10)  and  (12): 


Xj  > 


hk  +  b'jk 
2 


(23) 


However,  Constraints  (7)  and  (23)  eontradiet  eaeh  other,  implying  that  P'  is  not  lattiee  point  feasible, 
whieh  eontradiets  the  theorem’s  premise. 


4.  Equations  (13)  and  (15)  hold.  This  ease  is  identieal  to  Case  (3)  above. 

Thus,  if  n  77  =  0  and  we  obtain  a  semi-rounding  w'*  of  x'*  as  required.  This  eompletes  the  proof. 
□ 
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3.3  Main  Theorems 


We  now  arrive  at  the  key  result  of  this  paper. 

Theorem  3.2  Let  P  :  A  •  X  >  b  denote  a  G2SAT  polyhedron  and  x*  denote  a  half-integral  MFS.  If  P  is 
lattice  point  feasible,  then  it  contains  a  lattice  point  z  such  that  |  |z  —  x*  |  |oo  <  i-e.,  i  is  a  rounding  of'x.*. 

Proof:  We  prove  the  theorem  by  induetion  on  the  length  of  x. 

Base  Case:  Let  x  =  x  E  K.  If  x*  is  a  MFS,  there  exists  a  eonstraint  x  >  b  that  holds  with  equality  for 
X*.  Thus,  the  theorem  holds  trivially  for  z  =  x*. 

Induction  Step:  Let  us  assume  that  the  theorem  holds  for  all  veetors  x  of  length  up  to  n. 

Consider  the  ease  when  x  E  .  Sinee  P  has  a  MFS,  by  Theorem  (3.1),  it  has  one  with  half-integral 
entries.  Let  x*  =  •  •  •  ,  MFS  of  P.  If  x*  is  integral,  we  set  z  to  x*  and  we  are 

done.  So,  let  us  assume  that  x*  has  some  odd  half-integral  entries.  Note  that  if  two  variables  Xi  and  Xj 
appear  together  in  a  eonstraint  of  P  that  holds  with  equality,  either  both  x^  and  x*  are  integral  or  both  are 
odd  half-integral. 

Projeet  variable  Xn+i  out  of  P  using  Fourier-Motzkin  elimination  with  eoeffieient  normalization  (FM- 
CN).  Let  P'  :  A'  •  x'  >  b'  be  the  resulting  system,  where  x'  E  M” . 

Suppose  there  exists  a  lattiee  point  solution  y  =  {yi,y2,  ■  ■  ■  ,  yn-\-i)  of  P.  Thus,  y'  =  (yi,  y2,  •  •  •  ,  yn) 
is  a  lattiee  point  solution  of  P'. 

Consider  x'*  =  •  5  ^n)-  show  that  there  exists  a  rounding  zl  =  (zi,Z2,  ■  ■  ■  ,  Zn)  of 

x'*  whieh  satisfies  P'.  We  eonsider  the  following  three  eases: 

x'*  is  in  the  interior  of  P',  i.e.,  none  of  the  eonstraints  in  A'  •  x'  >  b'  hold  with  equality.  By  Corol¬ 
lary  3.1,  any  rounding  of  x'*  yields  a  lattiee  point  solution  zl  of  P'. 

Suppose  that  x'*  is  a  solution  of  P'  that  satisfies  some  eonsfraints  with  equality.  Suppose  that  for 
some  (A",  b")  C  (A',  b').  A"  •  x'*  =  b",  and  the  remaining  eonstraints  are  striet,  i.e.,  not  satisfied 
with  equality.  Sinee  x'*  is  a  MFS  of  A"  •  x'  >  b",  by  the  induetion  hypothesis,  we  ean  eonelude  that 
there  exists  a  lattiee  point  rounding  zl  of  x'*,  sueh  that  zl  is  a  solution  of  A"  •  x'  >  b".  Sinee,  by 
Corollary  3.1,  any  rounding  of  x'*  satisfies  fhe  sfriet  eonsfraints,  z'  is  also  a  lattiee  point  solution  of 
P'. 

It  is  possible  that  after  eoeffieient  normalization,  x'*  does  not  satisfy  P'.  By  Proposition  3.1,  there 
exists  a  semi-rounding  w'*  of  x'*  that  satisfies  P'.  Thus,  either  Case  (1)  or  Case  (2)  applies  with  x'* 
replaeed  by  w'*,  and  we  ean  obtain  a  rounding  z'  of  w'*  that  is  a  lattiee  point  solution  of  P'.  Finally, 
note  that  a  rounding  of  w'*  is  also  a  rounding  of  x'*,  sinee  integral  eomponents  of  x'*  are  preserved 
in  w'*.  This  eompletes  Case  (3). 

Thus,  we  ean  obtain  a  lattiee  point  solution  zl  of  P'  that  is  a  rounding  of  x'*. 

Sinee  P  is  G2SAT,  and  P'  is  obtained  from  P  using  FM-CN,  a  lattiee  point  solution  of  P'  ean  be 
extended  to  one  of  P.  Thus,  there  exists  an  integral  Zn+i  sueh  that  z  =  {zi^Z2t  ■  ■  ,  ZmZn+i)  is  a  solution 
ofP. 

To  eomplete  the  proof,  we  show  that  there  exists  sueh  an  integral  Zn+i  that  is  moreover  a  rounding  of 
Sinee  x*  is  a  MFS  of  P,  there  exists  a  subset  of  eonstraints  (A,  b)  of  (A,  b)  that  hold  with  equality 
at  X*.  The  value  of  Xn+i  is  eonstrained  only  by  the  values  of  other  variables  Xj  sueh  that  there  exists  an 
equation  in  Ax  =  b  in  whieh  Xn+i  and  Xj  appear  together.  Let  J  be  the  index  set  of  all  sueh  variables  Xj. 
We  now  show  that  there  exists  a  rounding  Zn+i  of  that  satisfies  Pi  :  Ax  >  b.  There  are  two  eases: 


Case  1: 

Case  2: 

Case  3: 
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1.  If  is  integral,  so  is  x*  for  all  j  E  J.  Thus,  Zn+i  =  satisfies  Pi,  and  we  are  done. 

2.  If  is  odd  half-integral,  so  is  x*j  for  all  j  E  J.  In  this  ease,  we  elaim  that  there  exists  a  eonsistent 

way  to  round  either  up  or  down,  so  that  the  result  satisfies  Pi.  Suppose  not,  i.e.,  there  exists 

eonstraints  that  foree  x’^j^i  to  be  rounded  up  as  well  as  down.  There  are  four  instanees  in  whieh  this 
might  oeeur: 

(a)  There  exist  eonstraints  Xn+i  —  Xi  >  b  and  Xj  —  Xn+i  >  6'  in  P  that  hold  with  equality  at  x*; 

furthermore,  zj  =  [x*]  and  Zi  =  Thus,  we  have  =  b+V ,  hut  Zj—Zi  <  b+V.  Sinee 

Xj  —  Xi  >  b  +  b'  is  a  valid  inequality  for  P,  this  means  that  z  does  not  lie  in  P,  a  eontradietion. 

(b)  There  exist  eonstraints  —Xn+i  —  Xi  >  b  and  Xj  +  Xn+i  >  6'  in  P  that  hold  with  equality  at  x*; 
furthermore,  Zj  =  \_x*j\  and  =  Kl-  This  ease  is  identieal  to  Case  (2a)  above. 

(e)  There  exist  eonstraints  Xj  —  Xn+i  >  b  and  Xj  +  Xn+i  >  6'  in  P  that  hold  with  equality  at  x*, 
with  Zj  =  \_x*j\ ■  Thus,  2a;^  =  6  +  6'.  Sinee,  x*^  is  odd  half-integral,  6+  6'  must  be  an  odd  integer. 
Moreover,  2zj  <  6  +  6'.  However,  sinee  2xj  >  6  +  6'  is  a  valid  inequality  for  P,  this  means  that 
z  does  not  lie  in  P,  a  eontradietion. 

(d)  There  exist  eonstraints  Xn+i  —  Xi  >  b  and  —Xn+i  —  Xj  >  6'  in  P  that  hold  with  equality  at  x*, 
with  This  ease  is  identieal  to  Case  (2e)  above. 

Thus,  there  exists  a  eonsistent  way  to  round  either  up  or  down  and  satisfy  every  eonstraint  in 
Pi.  Let  Zn+i  be  this  rounding. 

Applying  Corollary  3.1,  any  rounding  of  x*  satisfies  fhe  eonstraints  in  (A,  b)  \  (A,  b). 

Thus,  we  ean  obtain  a  rounding  z  of  x*  that  is  a  lattiee  point  solution  of  P. 

□ 

From  Theorem  (3.1)  and  Theorem  (3.2),  we  ean  eonelude  the  following  theorem. 

Theorem  3.3  Let  P  :  A  •  x  >  b  denote  a  G2SAT polyhedron  with  m  constraints  and  n  variables.  Then,  P 
has  enumeration  bound  n  ■  6niax- 

The  above  result  is  easily  generalized  for  arbitrary  G2SAT  formulas. 

Theorem  3.4  Let  (j)  denote  a  G2SAT formula  with  m  constraints,  n  variables,  and  let  6niax  be  the  maximum 
over  the  absolute  values  of  constant  terms  appearing  in  f.  Then,  f  has  enumeration  bound  n  ■  {bmax  +  1)- 

Proof:  If  (j)  has  a  satisfying  integer  solution,  that  solution  must  satisfy  one  of  the  terms  in  the  disjune- 
tive  normal  form  (DNF)  of  f.  Eaeh  term  in  the  DNF  representation  of  ^  is  a  G2SAT  polyhedron  in 
whieh  the  eonstant  term  in  any  eonstraint  has  absolute  value  at  most  6niax  +  1  (we  use  6niax  +  1  in  plaee 
of  6niax  to  aeeount  for  eliminating  negations  on  eonstraints).  It  follows  that  there  is  a  solution  to  f  in 

[-n  ■  (6max  +  1),  n-  (6  max  +  !)]■□ 


3.4  Approximation  Results  for  Optimization 

Consider  the  problem  of  optimizing  an  arbitrary  linear  funetion  over  a  G2SAT  polyhedron  P.  This  problem 
is  NP-hard  (minimum  vertex  eover  is  a  speeial  ease).  As  a  eorollary  of  Theorem  (3.2),  we  obtain  the 
following  theorem  showing  that  one  ean  approximate  the  optimal  value  to  within  an  additive  faetor. 
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Theorem  3.5  Let  P  =  {x  :  A  •  x  >  b}  denote  a  G2SAT  polyhedron  that  contains  a  lattice  point.  Let  the 
integer  linear  program  be  max{c  •  x  :  x  E  P}. 

If  the  optimum  value  is  finite,  solving  the  LP-relaxation  and  rounding  the  solution  can  yield  a  feasible 

Zl”-  |c'| 

lattice  point  that  approximates  the  optimum  to  within  an  additive  factor  of  A — ^ If  the  LP-relaxation 
is  unbounded,  so  is  the  integer  program. 


Proof:  If  the  optimum  value  v*  of  the  LP-relaxation  is  finite,  it  is  attained  at  a  MFS  x*.  Since  P  is  lattice 
point  feasible,  by  Theorem  3.2,  there  exists  a  lattice  point  z  in  P  such  that  such  that  ||z  —  x*||oo  <  5.  It 


->  ->  Z”-  k'l 

follows  that  c  •  z  is  within  ±  — —  of  v*,  and  hence  of  the  integer  optimum. 

If  the  LP-relaxation  is  unbounded,  so  must  the  integer  program,  since  P  is  lattice  point  feasible  [17].  □ 


Moreover,  an  approximate  solution  can  be  obtained  in  polynomial  time  in  the  following  three  steps: 

1.  Check  whether  P  is  lattice  point  feasible  using  Fourier-Motzkin  elimination  with  coefficient  normaliza¬ 
tion.  If  P  is  lattice  point  infeasible,  stop. 

2.  If  P  is  lattice  point  feasible,  solve  its  LP-relaxation.  If  it  is  unbounded,  we  conclude  that  the  original  IP 
is  also  unbounded.  Otherwise,  the  optimum  is  attained  at  a  MFS  x* . 

_  Ici 

3.  Round  X*  to  obtain  an  integer  solution  that  is  within  ±  — — —  of  the  optimum.  The  rounding  is  per¬ 
formed  as  follows.  For  each  variable  Xi  that  has  an  odd  half-integral  value  x*^,  we  check  whether  adding 
the  constraint  Xi  =  [x*]  to  P  preserves  lattice  point  feasibility.  If  not,  we  set  Xi  to  and  iterate, 
picking  another  variable  to  round,  until  we  have  obtained  a  feasible  integer  solution. 


It  is  easy  to  see  that  each  step  can  be  performed  in  polynomial  time.  Notice  that  if  lattice  point  feasibility 
is  preserved  by  setting  Xi  either  to  or  to  [a:|J ,  the  direction  of  rounding  can  be  chosen  heuristically  to 
obtain  a  tighter  approximation. 

Our  approximation  theorem  is  general,  in  that  it  applies  to  any  generalized  2SAT  integer  program,  in¬ 
cluding  non  0-1  programs  with  arbitrary  coefficients  in  the  objective  function.  However,  the  approximation 
factor  is  additive,  and  the  result  is  more  likely  to  be  useful  for  non  0-1  programs.  In  contrast,  the  results  of 
Hochbaum  et  al.  [13]  guarantee  a  2-approximation  for  G2SAT  integer  programs  expressed  as  a  minimization 
problem  where  the  objective  function  is  required  to  have  non-negative  coefficients. 


4  Experimental  Evaluation 

We  now  present  experimental  results  demonstrating  that  a  decision  procedure  based  on  the  solution  bound 
derived  in  this  paper  can  outperform  other  state-of-the-art  procedures. 

4.1  Implementation 

We  implemented  a  decision  procedure  that  operates  in  three  steps.  First,  given  a  G2SAT  formula  f,  it 
computes  the  enumeration  bound  n  •  (6max  + 1)  ■  Second,  it  translates  the  input  G2SAT  formula  to  a  Boolean 
formula  by  replacing  each  integer  variable  by  a  finite-precision,  signed  bit-vector  that  can  take  any  value  in 
the  range  [— n •  (6max  +  1)5  ti-  (6max  + 1)] ■  Arithmetic  and  relational  operators  are  then  encoded  as  arithmetic 
circuits  and  comparators.  Let  fibooi  denote  the  resulting  Boolean  formula.  Clearly,  fibooi  is  satisfiable  if  and 
only  if  (j)  is  satisfiable.  Thus,  fhe  final  step  consisfs  of  invoking  a  Boolean  safisfiabilily  (SAT)  solver  on 
fibooi-  Notice  fhaf  fhe  franslafion  fo  SAT  lakes  polynomial  lime  and  lhal  fhe  size  of  fbooi  is  polynomial  in 
lhal  of  <p. 
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The  main  reason  for  using  a  translation  to  SAT,  as  opposed  to  a  non-SAT-based  proeedure,  is  that  our 
benehmarks  possess  a  non-trivial  Boolean  strueture.  Also,  by  this  approaeh,  we  ean  leverage  the  reeent  ad- 
vanees  in  SAT  solving  (e.g.,  [16,  9]).  For  our  experiments,  we  employed  the  zChaff  satisfiability  solver  [16]; 
however,  note  that  any  alternative  SAT  solver  ean  be  employed  just  as  easily. 

4.2  Setup 

A  set  of  randomly  generated  G2SAT  formulas  was  used  for  the  experimental  evaluation.  A  G2SAT  formula 
ean  be  viewed  as  a  Boolean  eireuit  where  the  inputs  to  the  eireuit  are  G2SAT  eonstraints  rather  than  being 
Boolean  variables.  Eaeh  formula  was  generated  based  on  3  parameters:  the  maximum  number  of  variables, 
an  upper  bound  on  the  size  of  the  eonstant  term,  and  the  maximum  depth  of  the  eireuit.  We  varied  the 
maximum  number  of  variables  over  the  set  {40, 80, 160,  320, 640},  the  eonstant  term  upper  bound  over  the 
set  {16,  256, 4096, 65536, 1048576},  and  the  maximum  eireuit  depth  over  {6,  7,  8, 9, 10}.  For  eaeh  ehoiee 
of  these  three  parameters,  we  generated  a  formula  using  one  of  three  different  random  seeds;  the  seed  was 
used  in  generating,  at  eaeh  level  in  the  eireuit,  either  a  randomly  ehosen  Boolean  operator  or  a  G2SAT 
eonstraint.  The  variables  and  eonstant  term  in  eaeh  G2SAT  eonstraint  were  randomly  generated  as  well. 
Finally,  the  resulting  G2S  AT  formula  was  eonjoined  with  a  set  of  upper  and  lower  bound  eonstraints  on  eaeh 
variable,  where  the  bounds  were  randomly  seleeted  to  be  between  0  and  the  upper  bound  on  the  eonstant 
term.  This  last  operation  was  performed  in  order  to  generate  a  mix  of  both  satisfiable  and  unsatisfiable 
formulas.  Thus,  in  total,  the  benehmark  suite  eomprises  375  formulas,  of  whieh  202  are  unsatisfiable. 

We  eompared  our  proeedure  againsf  fwo  ofher  deeision  proeedures.  Bofh  are  based  on  a  eombinafion 
of  a  SAT  solver  wifh  a  solver  for  a  sysfem  of  infeger  linear  eonsfrainfs.  The  firsf  is  a  publiely  available 
fheorem  prover  ealled  CVC-Life  [7]  (fhe  version  available  as  of  Deeember  2004).  CVC-Life  uses  a  SAT 
solver  for  finding  Boolean  assignmenfs  fo  fhe  formula,  freafing  G2SAT  eonsfrainfs  as  Boolean  liferals.  For 
every  sueh  assignmenf,  if  deeides  fhe  feasibilify  of  fhe  eorresponding  eonjunefion  of  G2SAT  eonsfrainfs  by 
using  fhe  FM-CN  proeedure  (if  aefually  uses  fhe  Omega  fesf  [19],  whieh  speeializes  fo  FM-CN  for  G2SAT 
eonsfrainfs).  Defails  abouf  CVC-Life’s  operafion  ean  be  found  in  fhe  papers  by  Barreff  ef  al.  and  Ganesh 
ef  al.  [2,  3].  The  SAT  solver  used  by  CVC-Life  is  a  modified  version  of  fhe  zChaff  solver  used  by  our 
proeedure.  The  seeond  deeision  proeedure,  wriffen  by  Daniel  Kroening  (eurrenfly  af  ETH  Ziirieh),  works  on 
similar  prineiples  fo  CVC-Life,  exeepf  fhaf  if  uses  fhe  CPLEX  eommereial  opfimizafion  soffware  [6]  (version 
9.0)  insfead  of  fhe  EM-CN  proeedure.  This  proeedure  also  uses  fhe  zChaff  solver  as  ifs  SAT  solving  engine. 

Experimenfs  were  run  on  a  Linux  worksfafion  wifh  a  2  GHz  Penfium  4  proeessor  and  1  GB  of  RAM.  Our 
deeision  proeedure,  ealled  UCLID,  is  wriffen  mosfly  in  Moseow  ML,  a  dialeef  of  Sfandard  ML.  A  fimeouf 
of  600  seeonds  was  imposed  on  eaeh  run. 

4.3  Comparison 

Ligures  1  and  2  eompare  UCLID’s  fofal  time  (fime  for  bofh  eneoding  and  SAT  solving)  fo  fhaf  faken  by 
CVC-Life  and  fhe  CPLEX-based  solver  respeefively.  In  eaeh  plof,  fhe  y-eoordinafe  of  a  poinf  is  fhe  fime 
faken  by  UCLID,  and  fhe  x-eoordinafe  is  fhe  fime  faken  by  fhe  deeision  proeedure  we  eompare  if  againsf. 
UCLID’s  fofal  fime  is  dominafed  by  fhe  SAT  solving  fime.  Nofe  fhaf  fhe  X  and  Y  axes  are  on  differenf 
scales.  This  is  because  UCLID  finishes  wifhin  30  seconds  on  all  benchmarks  whereas  fhe  run-fimes  for  fhe 
ofher  solvers  are  spread  ouf  over  fhe  enfire  range  [0, 600].^ 

Lirsf,  consider  fhe  comparison  wifh  CVC-Life.  We  observe  from  Ligure  1  fhaf  CVC-Life  performs  worse 
fhan  UCLID  overall,  timing  ouf  on  95  of  fhe  375  benchmarks.  However,  nofe  fhaf  fhere  are  171  benchmarks 

^Details  of  experimental  results  described  in  this  section  are  available  at  http :  /  /www .  cs  .  emu  .edu/~uclid/2sat. 
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Figure  1:  Experimental  comparison  of  UCLID  versus  CVC-Lite  for  G2SAT  formulas.  Note  that  the 
seale  on  the  Y-axis  is  about  20  times  that  of  the  X-axis. 

on  whieh  CVC-Lite  outperforms  UCLID.  UCLID  eompletes  within  15  seeonds  on  all  of  these  benehmarks, 
and  within  5  seeonds  on  all  but  22  of  them. 

The  eomparison  with  the  CPLEX-based  solver  yields  similar  results,  as  one  ean  observe  in  Figure  2.  In 
faet,  the  CPLEX-based  solver  even  performs  worse  than  CVC-Eite,  timing  out  on  246  of  the  375  beneh¬ 
marks.  UCEID  is  outperformed  on  only  16  benehmarks,  on  all  of  whieh  it  terminates  within  7  seeonds. 


Eigure  2:  Experimental  comparison  of  UCLID  versus  CPLEX-based  solver  for  G2SAT  formulas.  Note 
that  the  seale  on  the  Y-axis  is  about  20  times  that  of  the  X-axis. 

We  further  analyzed  our  results  by  dividing  the  benehmarks  into  4  eategories,  with  eaeh  eategory  eom- 
prising  benehmarks  on  whieh  UCLID’s  time  falls  within  a  eertain  range.  Eor  eaeh  eategory,  we  eomputed 
the  pereentage  of  benehmarks  on  whieh  UCLID  outperforms  the  other  two  solvers.  This  data  is  displayed 
in  Table  1.  We  note  that  the  benehmarks  on  whieh  UCLID  is  outperformed  are  those  on  whieh  both  it  and 
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the  competing  solver  finish  within  a  few  seconds.  Note  also  that  UCLID  finishes  wifhin  5  seconds  on  over 
80%  of  fhe  benchmarks. 


UCEID  time  range 
(time  in  seconds) 

Number  of 

benchmarks 

%  of  benchmarks  on  which  UCEID  runs  faster 

CVC-Eite  prover 

CPEEX-based  solver 

[0,  5] 

315 

52.38 

95.24 

(5,  10] 

42 

54.76 

97.62 

(10,  20] 

15 

80.00 

100.00 

(20,  30) 

3 

100.00 

100.00 

Table  1 :  Comparing  UCLID  with  other  solvers  using  a  time- wise  break-up  of  benchmarks.  The  second 
column  indicates  the  number  of  benchmarks  on  which  UCLID’s  run-time  is  within  the  indicated  range. 


Thus,  we  can  conclude  that  the  enumerative  approach  presented  in  this  paper  can  greatly  outperform  a 
more  traditional  approach  based  on  combining  a  SAT  solver  with  a  constraint  solver.  The  main  reason  for 
this  seems  to  be  that  solvers  based  on  the  latter  approach  enumerate  several  SAT  assignments  that,  while 
satisfying  the  Boolean  skeleton  of  the  formula,  correspond  to  infeasible  systems  of  G2SAT  constraints.  On 
the  other  hand,  UCLID’s  encoding  adds  in  all  the  “G2SAT  information”  necessary  for  the  SAT  solver  to 
significantly  prune  its  search  space. 

5  Conclusion 

We  have  proposed  a  new  approach  to  deciding  the  satisfiability  of  Boolean  combinations  of  generalized 
2SAT  constraints.  The  central  insight  is  that  it  is  sufficient  to  search  for  bounded  solutions,  where  each 
variable  is  restricted  within  the  finite  range  [— n  •  (6max  +  1),  «  •  (&max  +  1)]-  The  solution  bound  we 
derive  improves  over  previous  results  by  an  exponential  factor.  The  key  step  in  our  derivation  is  a  novel 
result  for  G2SAT  polyhedra  on  finding  integer  solutions  by  rounding  minimal  face  solutions.  Experiments 
demonstrate  the  efficacy  of  a  SAT-based  decision  procedure  based  on  our  theoretical  results. 

It  would  be  interesting  to  extend  our  results  to  Boolean  combinations  of  linear  constraints  that  comprise 
mostly  of  G2SAT  constraints.  Previous  work  [21]  has  shown  that,  for  Boolean  combinations  of  mostly  dif¬ 
ference  constraints,  the  exponential  term  in  the  solution  bound  depends  only  on  the  number  and  coefficients 
of  the  non-difference  constraints.  It  is  still  open  as  to  whether  a  similar  result  can  also  be  obtained  for 
formulas  of  mostly  G2SAT  constraints. 
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